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Abstract 

In this study, approximation formulas for evaluating Landau constants are elab- 
orated by using the Boubaker Polynomials Expansion Scheme (BPES). Results are 
compared to some referred studies. 
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1 Introduction 

Landau's constants are defined (U-lll for all positive integers n, by: 

These constants were defined in relation with a set F of complex analytic functions / 
defined on an open region containing the closure of the unit disk D (D = {z : \z\ < 1}) 
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satisfying the conditions: 



d 



/(z)| z=0 = and -^/(z)| z=0 = 1. 



(1.2) 



If £(f) is the supremum of all numbers such that /(D) contains a disk of radius 1, and 
that / verifies the additional conditions: 



|/(z)| < l(|z| < 1) and f(z) = J^a k z k 

k=o 



then: 



L = inf{£(/):/GF} 



k=Q 



(1.3) 



(1.4) 



In this paper, a convergent protocol is proposed in order to give analytical expressions 
to the Landau's constants. 



2 Resolution process 

In concordance of the approximations performed by Watson Q, Zhao (H and Popa 1(3]], 
a sequence {co n }n^o is defined as: 

11 B 
cv n = G n - - ln(n + A) (7 + In 16) + —— ^, (2.1) 

7T 7T M + C 

where 7 is Euler's constant and A, B and C are unknown constants. 

The resolution process aims to find accurate approximation to the values of the pa- 
rameters A,B and C such that {co n } n ^o is the fastest sequence which would converge to 
zero. 



3 Approximation using The Boubaker Polynomials Expan- 
sion Scheme (BPES) 

3.1 Presentation 

The Boubaker Polynomials Expansion Scheme BPES [5 J- [24 J is a resolution protocol 
which has been successfully applied to several applied-physics and mathematics prob- 
lems. The BPES protocol ensures the validity of the related boundary conditions regard- 
less main equation features. The BPES is mainly based on Boubaker polynomials first 
derivatives properties: 
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and 



with 



N N 



k=\ 



(3.1) 



£ dB ik (x) 
k=l 



N dB± k (x) , 



0> E 

k=l 



N 

^=1 



(3.2) 



B 4(n+l)( r n) 

Several solution have been proposed through the BPES in many fields such as numerical 
analysis |0-|l8l, theoretical physics BH-IHl, mathematical algorithms fV2\ , heat transfer 
fl3f, homodynamic |[T4l[T5l , material characterization ||l~6Tl , fuzzy systems modeling |fT7fl- 
E2I and biology 03131 • 



3.2 Application 

The resolution protocol is based on setting A, B and C as estimators to the constants A, B 
and C, respectively: 



(3.3) 



where B^ are the 4/c-order Boubaker polynomials IH5B -[22 | , r k are B^ minimal positive 
roots, No is a prefixed integer, and dk\k=i,...,N are unknown pondering real coefficients. 

As a first step, the coefficients, the coefficients ^\k=i,...,N are determined through 
Falaleev approximation [yy]: 



A = 


l 

2N 




B = 


1 

2N 


E£lff B 4Jk(»lt)/ 


C = 


1 

2N 


E£ 1 $B 4fc (xr fc ), 



1 

7T 



ln|n + ^)+7 + lnl6 



(3.4) 



The BPES solution for A is obtained by determining the non-null set of coefficients 



£fc|fc=i,...,M) that minimizes the absolute difference Ajv : 



1 N ° - 5 

2N fc=i 



7T 



(7 + In 16) 



(3.5) 
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with 



3 f 1 N ° 
A k = 7 /. E *B ik (xr k )dx. 

k=l 



4 Jo 



Values of B and C are consecutively deduced from coefficients £?\k=i,...,N an d ^u\k=l,...,No 



which minimize the absolute difference A' Nq : 



a; 



Nn 



1 A/o i i No / -I No 



E x fc - -(t+ ini6) + ^ e #n ( 1 

fe=i 



7T ~ ' 2N — a 



with 



3 /-l N 1 N 

X k = 7~ E x%4k(xrk)dx and Y k = Z k = V B ik (xr k )dx. 
An Jo fc ti fc ti 



Hence the final solution is: 



B 



G M = iln(n + A) + i( 7 + lnl6)- 

7T 7T n + L 



with 



A 


= A = 


l 

2N 


v-N 

Mfc=l 


ZkB 4 k(xr k ), 


B 


= B = 


1 

2N 


V-N 


Zp±k{xr k ), 


C 


= C = 


1 

2N 


^N 


^B 4k (xr k ). 



4 Results, plots and discussion 



(3.6) 



(3.7) 



Numerical solutions obtained by the given method are gathered in Table CD along with 
precedent refereed approximations by Falaleev [25J and Brutman ||2~6ll . 

Plots of the BPES solution are presented in Fig. [TJ along with referred solutions 11251 
H. 

Figure |2] displays the obtained values along with those recorded by Falaleev [25j 
and Brutman ||26| . For accuracy purposes, error analysis has been carried out for the 
two referred datasets. Examination of the quadratic error plots (Fig. |2j) shows that the 
amplitudes of the error are more exaggerated according to Falaleev approximation [25J, 
particularly for low values of n. 

Moreover, Figure |2] monitors an obvious logarithmic profile in the range with a 
quadrature error which doesn't exceed 3.5% on the whole range (Fig. |2]). The obtained 
profile is in good agreement with the values recorded elsewhere [25J-[35J. 
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Table 1: Landau constants values. 



n 


Landau constants 


Quadratic Error vs. Exact 


Falaleev Approximation ] 25 1 


Brutman Approximation ] 26] 


BPES Approximation 


Ref. |25[ 


Ref. 1 26 1 





0.97463516 


0.97469795 


0.97473411 


9.7924E-9 


1.30765E-9 


1 


1.24433844 


1.24440124 


1.24449240 


2.37026E-8 


8.31042E-9 


2 


1.38820978 


1.38827257 


1.38841873 


4.36628E-8 


2.13632E-8 


3 


1.48693516 


1.48699795 


1.48719911 


6.96731E-8 


4.04659E-8 


4 


1.56218004 


1.56224284 


1.56249900 


1.01733E-7 


6.56187E-8 


5 


1.62299481 


1.62305761 


1.62336877 


1.39843E-7 


9.68215E-8 


6 


1.67403346 


1.67409626 


1.67446242 


1.84004E-7 


1.34074E-7 


7 


1.71800808 


1.71807088 


1.71849204 


2.34214E-7 


1.77377E-7 


8 


1.75663844 


1.75670124 


1.75717740 


2.90474E-7 


2.2673E-7 


9 


1.79108390 


1.79114669 


1.79167785 


3.52784E-7 


2.82133E-7 


10 


1.82216319 


1.82222598 


1.82281214 


4.21145E-7 


3.43585E-7 


11 


1.85047605 


1.85053884 


1.85118001 


4.95555E-7 


4.11088E-7 


12 


1.87647497 


1.87653777 


1.87723393 


5.76015E-7 


4.84641E-7 


13 


1.90050977 


1.90057257 


1.90132373 


6.62525E-7 


5.64244E-7 


14 


1.92285648 


1.92291928 


1.92372544 


7.55085E-7 


6.49896E-7 


15 


1.94373675 


1.94379954 


1.94466070 


8.53696E-7 


7.41599E-7 


16 


1.96333123 


1.96339403 


1.96431019 


9.58356E-7 


8.39352E-7 


17 


1.98178915 


1.98185195 


1.98282311 


1.06907E-6 


9.43155E-7 


18 


1.99923515 


1.99929795 


2.00032411 


1.18583E-6 


1.05301E-6 


19 


2.01577445 


2.01583725 


2.01691841 


1.30864E-6 


1.16891E-6 
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Figure 1: Solution plots 
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Figure 2: Absolute error plots 

5 Conclusion 

In this study, approximation formulas for evaluating Landau constants have been pre- 
sented and discussed. The proposed estimate has been compared with two other es- 
timates which are of special importance in approximation theory. Results have been 
favorable for the performed method in terms of both convergence and accuracy. 
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